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ABSTRACT 

Confexf. In the last years several exoplanets have been discovered that orbit one component of a compact binary system (separation 
< 50 astronomical units), the probably best-known case is y-Cephei. So far, all attempts to explain the in-situ formation of these planets 
has been unsuccessful, in part because of the strong gravitational perturbations of the secondary star on any initial planetesimal swarm. 
Aims. Here we test whether planetary bodies in compact binaries, in particular y-Cephei, could have originated from a close encounter 
with a passing star, assuming initial configurations for the stellar system suitable for planetary formation. In other words, we analyze 
whether the orbital configuration of the current binary system might have been generated after the formation of the planet, and as a 
consequence of a close encounter with a third star in hyperbolic orbit. 

Methods. We performed a series of time-reverse N-body simulations of stellar scattering events in which the present-day configuration 
of y-Cephei was used as the initial condition plus a hypothetical third star as an impactor. We analyzed which configurations and 
system parameters could have given birth to the current system. 

Results. Depending on the maximum impact velocity allowed for accretional collisions, we find that between ~ 1% and ~ 5% of 
stellar encounters correspond to an "original" system in which planetary formation around the primary star is not inhibited by the 
secondary, but is acceptable within the classical core-accretion scenario. Thus, although not highly probable, it is plausible that stellar 
encounters may have played a significant role in shaping these types of exoplanetary systems. 

Key words, planets and satellites: formation - stars: binaries: y-Cephei. 



1. Introduction 

Studies of planetary formation in binary stellar systems are an 
intriguing and challenging task. Although the majority of stars 
are members of such systems (Duquennoy & Mayor 1991), only 
a handful of exoplanets have so far been detected around com- 
pact binaries (i.e. stellar separation below ~ 50 AU). Even so, 
their very existence appears to be inconsistent with models and 
theories that work for single stars, which pushes some crucial 
parameters to extreme values. 

The case of the y-Cephei binary is particularly interesting 
because it is one of the most extreme systems that harbors a giant 
planet. The stellar system is formed by a primary star with a 
mass of Ma = 1.6M Q and a secondary with a mass of Mb = 
O.34M . The Mb body orbits around Ma with a semimajor axis 
of ciab = 18.5 AU and an eccentricity of cab = 0.36 (Hatzes et 
al. 2003). 

The giant planet that orbits Ma has been the object of 
many studies over the past decade (e.g. Marzari & Scholl 2002; 
Thebault 2004; Thebault 2006; Haghighipour & Raymond 2007; 
Paardekooper et al. 2008; Beauge et al. 2010; etc.) and we have 
yet to find a suitable physical scenario to explain its existence. 
Artymowicz and Lubow (1994) showed that the circumplanetary 
disk is tidally truncated by its companion. This truncation oc- 
curs at a location close to the outer limit for dynamical stability 
of solid particles, and therefore safely lies beyond the position 
of all detected exoplanets. It nevertheless poses a problem, es- 
pecially for giant planet formation, because it deprives the disk 
from a large fraction of its mass. Another consequence of disk 
truncation is that it shortens the lifetime of the disk, and in turn 
the timespan for gaseous planet formation (Cieza et al. 2009). 
To complicate the scenario even more, Nelson (2000) suggested 



that for an equal-mass binary of separation 50 AU, temperatures 
in the disk might stay too high to allow grains to condense. 

Even if grains can condense, binary perturbations might im- 
pede their growth by increasing the impact velocities beyond the 
disruption limit (e.g. Marzari & Scholl 2000). Although gas drag 
helps to introduce a phase alignment in the pericenter and, thus, 
to decrease the relative velocities of nearby bodies, this appears 
to work only for equal-size planetesimals (e.g. Thebault et al. 
2006) and circular gas disks. Eccentric and even precessing gas 
disks complicate the picture still more (Paardekooper et al. 2009; 
Marzari et al. 2009; Marzari et al. 2012), although it may indeed 
work in our favor in the outer parts of the disk (Beauge et al. 
2010). 

In recent years several studies have been performed on the 
role of mutual inclinations between the planetesimals, the gas 
disk, and the binary companion (Marzari et al. 2009b; Xie et 
al. 2010; Thebault et al. 2010; Midler & Kley 2012; Zhao et al. 
2012.). However, even this additional degree of freedom does not 
appear to be a solution and, in the particular case of y-Cephei, 
does not seem to favor planetary formation for bodies at a > 2 
AU, which is expected to be the breading ground for the planet 
detected in this system. 

In a recent review of the state of the art in this problem, 
Thebault (2011) wondered whether the whole idea of core- 
instability planetary formation is feasible in such systems, or 
if we must resort to gas-instability models (Duchene 2010). 
However, as Thebault (201 1) points out, "the disk instability sce- 
nario also encounters severe difficulties in the context of close 
binaries, and it is too early to know if it can be a viable alterna- 
tive formation channel". 

So, the question remains: How did the planet manage to form 
around the star y-Cephei-A? Here we revisit an alternative route, 
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based on the concept of stellar scattering. We assume that planet 
m p formed (by one of the accepted models of planetary forma- 
tion) around Ma, but in a different dynamical configuration in 
which the star Ma is single or the companion Mb is sufficiently 
wide to allow the formation of the planetary system. 

The scenario explained is not new. Pfahl (2005) and 
Portegies Zwart & McMillan (2005) discussed such a mecha- 
nism to explain the origin of a planet then believed to exist in 
the HD 188753 system (Konacki 2005), although the existence 
of this body was subsequently questioned and is so far uncon- 
firmed. 

A slightly different analysis was later presented by Marzari 
& Barbieri (2007a, 2007b), were the authors studied the survival 
of circumstellar planets in dynamically unstable triple stellar 
systems or in binaries that suffered a close encounter with a third 
(background) star. These papers explored, for the first time the 
possibility that the current configuration of the y-Cephei stellar 
components may not be primordial, and that the exoplanet may 
have been formed under different circumstances. 

In this paper we again address the proposal of Marzari & 
Barbieri (2007a), and study the effects of hyperbolic fly-bys of 
background stars in a binary system similar to y-Cephei. The 
idea is to see under which initial conditions and system param- 
eters it is possible to transform a primordial accretion-friendly 
stellar system to the current (accretion-hostile) configuration. 
Since our is dedicated to a single binary system, we are able to 
explore a larger set of parameters and situations, as well as dif- 
ferent types of outcomes. We also include a statistical analysis 
of the results considering distribution probabilities for impacting 
velocities and stellar masses. 

The work is organized as follows. In Section 2 we review 
the idea of gravitational scattering, with special emphasis on hi- 
erarchical three-body systems. Section 3 presents our N-body 
code and some test runs to analyze the importance of the free 
parameters in the scattering outcomes. Our main numerical sim- 
ulations are discussed in Section 4, while a probability analysis 
of positive outcomes is presented in Section 5. Conclusions fi- 
nally close the paper in Section 6. 



2. Gravitational scattering in hierarchical 
three-body systems 

The gravitational interaction between a star and a binary sys- 
tem is a classical problem in stellar dynamics (e.g. Heggie 1975; 
Hills 1975; Hut & Bahcall 1983; Heggie & Hut 1993; Portegies 
Zwart et al. 1997; etc.), although it is focused on the relationship 
between binaries and the dynamical evolution of open clusters. 
The same principles, however, may also be applied to other prob- 
lems, such as the origin of exoplanets in multiple stellar systems 
(Pfahl 2005, Portegies Zwart & McMillan 2005) and the dynam- 
ical evolution of "cold" exoplanets with orbits very distant from 
the central star (e.g. Veras et al. 2009). 



2. 1 . The hyperbolic two-body problem 

We assume a compact stellar binary system with masses iua and 
rriB, whose orbital separation is small compared with the initial 
distance of an incoming star mc- In this scenario, and at least 
for the initial conditions of the hyperbolic encounter, we may 
approximate the two binary components by a unique body of 
mass rtiA + m B (see Hut & Bahcall 1983), and refer the orbit of 
the impactor with respect to the barycenter of the binary. 



In classical scattering problems, the outcome of the fly-by 
may be defined by two fundamental parameters: the infall veloc- 
ity of the incoming body at infinity Voo = |r(f — > -oo)| and the 
impact parameter b, defined as the minimum distance between 
the two bodies if there were no deflection. Alternatively, it is pos- 
sible to substitute b with the minimum distance r m ,„ = a(l - e) 
between the impactor and the center of mass of the pair niA+m B . 

Following the approach developed by Hills (1975), we can 
approximate the total orbital energy of the three-star system, 
prior to the encounter, by 



E, n , — Ear + Er — — 



K 1 2 

2a A B 2 



(1) 



where both k and // are two mass factors defined according to 



k = Q(m A + m B )m c 



V = 



(m A + m B ) m c 
iua + m B + trie 



(2) 



In equation (Q]) E A b is the binding energy of the binary and Ec is 
the kinetic energy of the impactor. oab is the semimajor axis of 
m B with respect to ma- We assume that initially mc is sufficiently 
far from the other two that we can neglect its potential energy. 
The same calculation can also be made after the scattering event 
has taken place, now in terms of the new semimajor axis of the 
binary (a' AB ) and the outgoing velocity of mc (which we will 
denote by v^). 

We can then use the conservation of E tot to write (see Heggie 
1975) 



~V v 'l, = ~A"'L + ae ab, 



(3) 



where AE A b is the change in the binding energy of the binary. 
The minimum velocity v c of the impactor necessary to disrupt 
the binary is such that AEab = E A b or, equivalently, E tot = 0. 
This gives (Hut & Bahcall 1983) 



k 1 QmAmsimA + m B + mc) 1 



Ha AB 



(m A + m B )m c 



cab 



(4) 



It will prove useful to write the incoming velocity of the im- 
pactor Voo in units of v c 



(5) 



which gives an adimensional expression for the incoming veloc- 
ity. Furthermore, v«, can be written in terms of the semimajor 
axis of the hyperbolic orbit as: 



Q(m A + m B + m c ) 



(6) 



where a is the the semimajor axis of mc with respect to the center 
of mass of the binary m,A +m B . Recalling expression ©, we may 
now write the relationship between a and v* as 



-ClAB' 



m c (m A + m B ) 1 



m A m B 



v* 1 



(7) 



Similarly, we can obtain the eccentricity of the hyperbolic orbit 
of mc as 



1 



(8) 



Finally, if we define 6o as the value of the eccentric anomaly 
at the beginning of the simulation (when r = ro), we can cal- 
culate the initial position and velocity of mc using the usual ex- 
pressions for hyperbolic motion (e.g. Butler 2005). 
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2.2. Complete set of variables for the scattering problem 

The set (v„ r min , eo) defines the initial conditions for the impactor 
m c , with respect to the center of mass of m A + m B , assuming 
that the differential gravitational interactions from both compo- 
nents do not significantly perturb the two-body hyperbolic orbit 
(see Hut & Bahcall 1983 for a similar approach). However, even 
if this condition is not satisfied, and the orbit of the incoming 
star is seriously affected by the two central bodies, the values of 
(v*, r min ) are still indicative of the "unperturbed" motion of m c , 
and therefore may still be referred to as variables of the problem. 

To complete the description of the system, we must give the 
initial conditions for each of the binary components at a given 
time. These can be uniquely specified by the set of m^-centric 
orbital elements of (a B , e B , I B , ^b, wg, SIb), where the incli- 
nation I B is given with respect to the orbital plane of the incom- 
ing impactor. Since the value of e is already a variable of the 
system, it appears less complicated to specify the value of the 
mean longitude A B at the time of maximum approach of the im- 
pacter (i.e. r = r m i„). From the two-body solution we can then 
calculate the value of the mean longitude at the initial time e = eo 
and calculate the position and velocity vectors of m A and m B at 
the beginning of the simulation. 

As with the initial conditions for mc, it is not expected 
that the orbits of the binary stars remain constant in time. 
Accordingly, the values for the orbital elements adopted as ini- 
tial conditions are also "unperturbed" values, in the same way 
as the impact parameter b is the unperturbed minimum distance 
between the projectile and the target. 

Last of all, and since our ultimate aim will be to study the 
exoplanet detected around y-Cephei-A, we will also assume that 
a small planetary mass m p orbits the central star m A . The initial 
m^-centric orbit of the planet will be circular and specified by 
values of (a p ,I p ,A p ,Cl p ), where A p is once again given at the 
time when r c = r min . 

To summarize, we need to set initial values for 16 indepen- 
dent variables in order to get a full description of the scattering 
event. These are 

(m A ,m B ,m c ,m p ) + (r min ,v*) + (9) 
(a B ,e B ,I B ,A B ,LO B ,n B ) + (a p ,I p ,A p ,Q. p ). 

We recall that the mA-centric orbital elements of m B are assumed 
to correspond to the point at which the impactor has its maxi- 
mum approach if the impactor had zero mass. 

Of course, previous knowledge of the dynamical system im- 
plies that some variables are already known, such as most of the 
masses and the main orbital elements. However, this still leaves 
us with a significant number of degrees of freedom to explore. 

3. The numerical code 

Once we set the values for the stellar masses and initial condi- 
tions, we integrated the equations of motion using an N-body 
code to simulate the scattering event. Dynamically, the problem 
may be characterized as a hierarchical (or nested) full four -body 
problem, in which all mutual gravitational interactions are taken 
into account. 

For the simulations described in this work we employed a 
code constructed around a Bulrisch-Stoer integrator with an er- 
ror tolerance of = 13. Since the scattering event is expected 
to be extremely sensitive to the initial conditions, especially if 
the impact parameter is small, we performed a series of prelim- 
inary runs that included back integrations from the final config- 



urations. In all cases there was no significant difference between 
the two branches of the same trajectory. 

Our working hypothesis is that the current y-Cephei system 
could be the outcome of a stellar scattering event that either 
modified the orbit of the secondary star or caused an exchange 
of massive bodies. To study whether this is possible, we em- 
ployes the concept of reverse integration. In other words, our 
initial conditions were the present day orbital configuration of 
the y-Cephei system plus an additional star mc in an arbitrary 
outward-bound hyperbolic orbit. We then integrated the system 
backwards in time to reproduce the hypothetical scattering event 
and find the "initial conditions" that led to the present system. It 
is important to stress that this procedure can only be employed in 
purely conservative systems where the equations of motion are 
time-reversible. Therefore, we did not consider the effects of any 
non-conservative force such as gas drag of tidal effects between 
the bodies. 

We considered that a simulation yields a positive result when 
the following two conditions are met: 

- The planet m p remains bounded to m A in an orbit with 
moderate to low eccentricity (e.g. e p < 0.1). 

- The region around m p becomes compatible with planetary 
formation according to standard core-accretion scenario. 
This could either be a consequence of a significant increase 
in the semimajor axis of the secondary binary component, or 
a stellar exchange in which m A becomes a single star. 

Since the equations of motion of the N-body problem are 
time-reversible, the final outcomes of our simulations would 
then constitute possible initial conditions that could explain the 
current planetary system. 

3.1. Possible outcomes 

Depending on the relative magnitudes between Voo and AE AB , 
we can list five possible outcomes of a hyperbolic fly-by: De- 
excitation, excitation, resonance, ionization, and exchange. The 
reader is referred to Heggie (1975) for more details, including 
explicit conditions for each type. 

Not all possible outcomes are relevant for the case at hand. 
For example, resonance is a transitory state, at least in the point- 
mass approximation, and following complex interactions be- 
tween all three stars, at least one will leave the system in a hy- 
perbolic orbit. Ionization implies that all three stars acquire un- 
bounded motion with respect to the other two; the original binary 
is then destroyed, but not replaced. Since in our case we worked 
with a reverse temporal axis, this would mean that the present- 
day y-Cephei is the outcome of a triple encounter of isolated 
stars, a highly improbable situation, to say the least. 

Since our aims are more restrictive and focused on finding 
possible formation mechanisms for y-Cephei-A-b, we prefer to 
redefine the possible outcomes according to what Pfahl (2005) 
calls "scenarios". These are: 

- Type I: The impactor m c captures the secondary component 
m B during the encounter, leaving behind m A as a single star 
with the planet still bounded in a quasi-circular orbit. 

- Type IT. While there is no exchange in the memberships, the 
binary m A + m B acquires a sufficiently wide orbit to allow 
planetary formation in the region around m p . 
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Fig. 1. Examples of scattering events that lead to all three types of outcomes. Black curves correspond to wa, red to mg, and light 
green to the impactor mc- 



- Type III: The incoming star mc changes place with m# where 
this star is expelled from the system. The new binary iha +mc 
is sufficiently wide to allow planetary formation around the 
location of m p . 

Both types I and III are cataloged as exchange in Heggie's 
nomenclature, while type II corresponds to an excitation. 
Marzari & Barbieri (2007a) only discussed type-II outcomes. 

In addition to the types listed before, we also found a few 
cases of three other types of outcomes. Although these were very 
rare and have only a negligible effect on the statistical analysis 
performed later on, we will nevertheless include them here for 
completeness and for a full understanding of the complexity of 
the problem. These will be referred to as types IV to VI, and 
correspond to the following scenarios: 

- Type IV: The impactor mc not only captures the secondary 
component m# during the encounter (which would be 
cataloged as type I), but also the planet that originally 
orbited m^. The final configuration then includes a new 
binary mg + mc with the planet orbiting mc. 

- Type V: The impactor mc carries away both m^ and its planet 
m p , but in the final configuration the planet now orbits mc- 

- Type VI: The passage of the impactor mc only captures the 
planet from the original binary, leaving behind a planet-less 
ma + mg system. 

As mentioned in the previous section (see equation (f2]i), we 
need to specify a total of 16 parameters/variables. Of these, how- 
ever, several will be fixed to represent the y-Cephei system. For 
the masses we adopted the values 



mA — 1.6to q ; mg — 0.4m Q ; m p — 1.7»?j up , 
while for the -centric orbit of ms we assumed: 
a B = 18.5 AU ; e B = 0.36. 



(10) 



(ID 



Finally, we assumed for the TO^-centric semimajor axis of the 
planet m p a p - 2.1 AU. The remaining ten parameters where 
considered to be unknown and free to modify in our series of 
runs. 



3.2. Test runs 

Not all variables are created equal, and some are expected to 
be more influential in determining the outcome of the scattering 



event. The aim of these preliminary series of runs is to estimate 
how each parameter affects the dynamics and what may be pref- 
erential values. 

We began with three series of runs, each with fixed values of 
the impactor mass. The values adopted for each run are mc = 
0.2 (series 1), mc = 0.4 (series 2) and mc = 1.0 (series 3), 
with values given in solar masses. All orbits are assumed to lie 
on the same plane (i.e. /, = 0, and Q, indeterminate) and the 
remaining angles are taken equal to zero. At the beginning of 
each integration, the impactor is set at a distance from the binary 
equal to ro = 5QQb, where Qb = Ob(1 + es) is the apocentric 
distance between both components. Since the orbital elements 
of the planet and binary are given at r — r m - ln , the outcome of the 
scattering event is not dependent on the value of ro, as long as it 
is chosen to be sufficiently large. 

Thus, in each series we have only two free parameters left: 
(fmin,v*). For each series we chose a total of 2.5 x 10 5 values 
chosen from an equispaced 500 x 500 grid within the intervals: 



1.1 < v-oAw < 10 
0.1 < r min /Q B <4, 



(12) 



where vo is the magnitude of the initial velocity at r = ro and v esc 
the escape velocity at that point. Thus, we consider initial speeds 
ranging from almost parabolic orbits to highly hyperbolic fiy- 
bys. After some simple algebraic calculations, we find that the 
corresponding interval in v* is approximately 0.08 < v* < 1.7. 

Figure [2] shows the post-scattering -centric orbital ele- 
ments of the secondary star of the resulting binary system. Since 
the central star becomes isolated in type-I interactions, only type 
II (left) and type III (right) are analyzed. In the first case, mj 
remains bounded to wa, so the eccentricity and semimajor axis 
plotted are those of m^. In type-Ill outcomes, however, the final 
binary is now ;«a + mc, so the orbital elements are those of mc 
with respect to niA- Black squares show those runs in which the 
planetary orbit remained almost circular, while those in which 
m p was severely perturbed are shown in brown. 

In most cases the resulting distribution is V-shaped, roughly 
delimited on the left by apocentric distances Q > 25 AU, which 
corresponds to the apocentric distance of the original binary sys- 
tem. The black squares, on the other hand, are mostly located 
with q > 12 AU, which is also the pericentric distance of the 
original system. It therefore appears that if the new binary has 
a lower value of q, roughly independent of the eccentricity, its 
gravitational perturbations on the planet would be sufficient to 
increase its eccentricity beyond the imposed limit e p - 0.1. Of 
course this analysis is qualitative and does not rigorously rep- 



J.G. Marti and C. Beauge: Stellar scattering and the origin of the planet around y-Cephei-A 



5 



Type II Type III 




10 100 1000 10 100 1000 




a c 



Fig. 2. Distribution of the final eccentricity and semimajor axis 
of the secondary binary with respect to the central star m^, after 
the scattering event. In the left-hand plots (type II) the secondary 
star is still mg, while in the right-hand plots it is mc- Brown 
squares show all runs, while those that preserved the planetary 
orbit with e p < 0. 1 are superimposed in black. The orange lines 
show level curves of typical impact velocities (see text for de- 
tails). 

resent every solution, but it does give an overall picture of the 
behavior. 

Not all black squares are categorized as positive solutions. 
As mentioned previously, the final system must also permit for 
planetary formation, i.e., allow for an in-situ accretion of m p . 
Since the present-day y-Cephei system is too compact to al- 
low for accretional collisions in the planetesimal swarm (e.g. 
Thebault et al. 2004; Paardekooper et al. 2008), we must ask 
about the minimum distance between binary components at 
which this may be avoided. 

Unfortunately, there is no easy answer to this question, but 
we may obtain a rough idea from previous numerical studies 
of planetesimal dynamics in binary systems. Thebault et al. 
(2006) presented a detailed calculation of the impact velocity 
of different-sized planetesimals located at a = 1 AU in generic 
binary systems with iua = 1 and mg = 0.5, both in units of solar 
mass. For every point in a grid of values of (as, eg) these au- 
thors calculated the average impact velocity Av between a 2.5 
km and a 5 km solid body, and drew level curves of constant Av. 
According to the authors, accretional collisions are guaranteed 
for Av < 10m/ s, denoting regions where the encounter speeds 
are practically unaffected by the secondary star. Collisions with 



10 < Av < lOOm/s are less certain, but may still allow for plan- 
etary formation even if some erosion is expected. 

Since the impact velocity at any given location is primar- 
ily dictated by the forced eccentricity ej excited by the sec- 
ondary star on the planetesimal swarm, we can use the results 
by Thebault et al. (2006) to extrapolate their findings to other 
binary systems. As determined by Heppenheimer (1978), a first- 
order expression (in the masses) for the forced eccentricity at a 
given semimajor axis a may be written as 



A second-order expression was calculated recently by Giuppone 
et al. (2011), which depends explicitly on the stellar masses. 
However, the difference between both values is only relevant for 
high values of e/, which is beside the purpose of the present 
study. 

Using expression CfT~3T > of Thebault's system we can calculate 
(for our case) the critical value of the forced eccentricity associ- 
ated to Av = 10m/ s and to Av = 100m/s. We can then use these 
as reference values and, once again applying ([T3T >. estimate the 
values of ag and eg that give the same value of e/ for a given 
semimajor axis a. Since the forced eccentricity increases with a, 
cib and <?b, in principle we should expect accretional collisions 
for lower values of all these parameters. 

To complete this setup, we must specify a value for a. Since 
the exoplanet detected in y-Cephei has a semimajor axis of 2.1 
AU, even an in-situ formation with little planetary migration re- 
quires a sufficiently large feeding zone that would extend be- 
yond the snow line to accrete the necessary mass. For the present 
study, we chose a — 4 AU; in other words, we required that 
the entire region up to 4 AU satisfies the conditions established 
in Thebault et al. (2006). Although this upper limit may appear 
rather arbitrary, we consider it as illustrative of the dynamics and 
not as an actual proposal for the formation process itself. 

The orange curves in each frame of Figure |2]show the values 
of the eccentricity and semimajor axis of the binary system that 
yield (within the approximation described above) the two values 
of the impact velocities mentioned before. All binaries located 
to the right and below each curve satisfy this ad-hoc condition 
of accretional collision. 

Figure [3] focuses on the type-II outcomes for ;«c = 0.4M Q . 
The top frame shows the distribution of orbital parameters of 
the binary as a function of the minimum distance r m ,„ of the 
impactor, in units of Qb- For illustration purposes, we have 
also drawn the curves of qs = «b(1 - <?/,) = 12 AU and 
Qb = ob(1 + <?/,) = 25 AU, both corresponding to the origi- 
nal values of the system. Thus, the original orbit is located at the 
intersection point between both curves. 

Since distant encounters (i.e. large values of r„„„) imply 
lower perturbations, the final configuration remains close to the 
original orbit. However, as the distance decreases, the resulting 
orbits spread significantly, most of them increasing both qs and 
Qb- 

A comparison with the curves of constant Av also shows that 
the positive solutions occur primarily for impactors whose min- 
imum distance lies between one and three times the apocentric 
distance of the original binary. Closer encounters are too destruc- 
tive, while more distant fly-bys are not able to modify the system 
sufficiently to overcome the planetary formation constraints. 

The lower plot shows the distribution of final orbits as a func- 
tion of the velocity at infinity of mc, in units of the escape ve- 
locity of the system. Low incoming speeds cause a more or less 
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Fig. 3. Distribution of the eccentricity and semimajor axis of m B 
with respect to m^, separated according to the values of r„,j„/QB 
(top) and v* (bottom). The orange curve on the right corresponds 
to q B = 30 AU. 
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Table 1. Percentage of positive outcomes per event type and im- 
pactor mass. 



random spread of solutions around the original configuration, 
mainly as a consequence of the resonance interactions. Higher 
speeds generate more ordered distributions. Once again, posi- 
tive solutions are mainly found for intermediate values, between 
two and six times the escape velocity. 

Finally, the percentage of positive results are summarized in 
Table Q] were we have extended the series to include a total of 



five values of rric between O.2M and 1M Q . For most impactor 
masses, type I positive solutions are the most common because 
they contain no restrictions to the orbit of the outgoing binary. 
We reecall that in this type of outcome the central star niA be- 
comes isolated and the new binary is composed of m B + mc. All 
other types of solutions are rare, as already noted by Marzari & 
Barbieri (2007a), usually less than 1% of all scattering events. If 
we consider the very restrictive Av < 10 m/s limit, only a few 
0.01% of the encounters yield final configurations suitable for 
planetary formation around niA- 



4. Full numerical simulations 

We now repeated the simulations, this time with random values 
of all angular variables {I B , A B , tj b , Q b , A p , m p ). In other words, 
we extended the experiments to the 3D case and considered arbi- 
trary orientations of the orbits of both the binary and the planet. 
However, in all cases we considered the planet initially in the 
same orbital plane as the secondary, i.e. I p = I B and £l p = Q B . 
One of the aspects we wished to study is the change in the mu- 
tual inclination between the planet and the binary perturber after 
the passage of the impactor. It was already noted in Marzari & 
Barbieri (2007ab) that high inclinations between the impactor 
and the binary may lead to greater disruptions of the system, so 
these configuration may cause large variations in I p , leaving be- 
hind planetary systems with high inclinations with respect to the 
plane of the binary. 

Figure |4] shows the eccentricity and semimajor axis distri- 
bution of the resulting binaries, where the left-hand plots show 
the type-II outcomes and the right-hand graphs correspond to 
type III. As with Figure [2] the symbols in brown show all the 
outcomes while those in black maintained the planet in an al- 
most circular orbit (e p < 0.1). Compared with Figure|2] we see a 
greater proportion of positive results for type-II outcomes, but a 
significant decrease in positive results for type III (see Table|2]i. 
In fact, the number of Type III outcomes with relative encounter 
speed Av < 10 m/s expected for small planetesimals is less than 
0.01% for all impactor masses. Results also show a serious re- 
duction in the positive outcomes of type-I encounters, which are 
now typically only half those found in our test runs. 

Finally, Figure [5] shows the distribution of final inclinations 
imut of the planetary orbit with respect to that of the binary com- 
ponent. In type-II outcomes (left-hand plots), the secondary re- 
mains the original companion m B , therefore i mut measures the 
perturbation of the fly-by on the original system. Although most 
orbits remain relatively aligned, especially those corresponding 
to positive results, we also note a significant amount of highly 
inclined orbits, some even retrograde with respect to the binary 
orbit. 

The right-hand plots in Figure [5] correspond to type-Ill out- 
comes, in which the new binary system is now composed of 
Ma + mc- Because the new companion is captured and not pri- 
mordial, in principle there should be little correlation between 
the orbits of m p and mc; consequently the values of i mut should 
be fairly random. However, the results of the simulations, es- 
pecially those classified as positive, show a marked correlation. 
This seems to indicate that only the fly-bys that occur near the 
orbital plane of the original binary lead to systems suitable for 
planetary formation. The difference in the plots between the 
black and red curves gives the impression that most of the en- 
counters with large I B lead to highly eccentric planets, and posi- 
tive results occur primarily with I B ~ 0. 
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Fig. 4. Same as Figure[2] but for random orientations and inclina- 
tions of the orbits. Top plots show results for an impactor mass of 
mc = 0.4M o , while the bottom plots correspond to mc = IM Q . 
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Table 2. Percentage of positive outcomes per event type and im- 
pactor mass for full 3D simulations. 



5. Distribution functions for the free parameters 

To gain a more accurate understanding of the probability to form 
the y-Cephei system through the process explained in these pa- 
per, we must weight the results showed in the previous sections 
by assigning a distribution function (DF) to each of the free pa- 
rameters, r m ,„, v* and mc- 

As mentioned earlier, the parameter r mi „ appears more suit- 
able than the impact parameter b for describing the two-body 
dynamics of the scattering problem. However, to find the corre- 
sponding DF, we must first relate r,„,„ with b and v* . 

For a hyperbolic keplerian orbit we may write 



r,nin = -a(e - 1) 
b = -a Ve 2 - 1 

a ' 



(14) 



where a and e are the semimajor axis and eccentricity of the orbit 
of the incoming body mc and A = ciab mc(mA + m^jirn/^B)- 



Introducing the last two expressions of (26) into the first, and 
after some simple algebraic manipulations, we obtain 



,{b,v*) = — 




(15) 



Following Hut & Bahcall (1983) we assumed a distribution 
of close encounters uniform in b 2 . To estimate the distribution 
probability of v*, we first need to determine the DF of the rel- 
ative encounter velocities between a single star and a binary. 
Again following Hut & Bahcall (1983), the thermal distribution 
function of velocities for stars of equal mass m is 



nj \kT) V 6XP \ 2kT V /' 



(16) 



where k is Boltzmann's constant and T is a temperature-like pa- 
rameter in analogy with a velocity DF of gas particles. Below 
we will express these variables in terms of the thermal velocity 
dispersion v>/,. 

Binaries formed from equal-mass stars have a combined 
mass 2m, and a thermal DF 



where T -T for equipartition of translational energy, and T 
27 if the binaries have the same velocity dispersion as single 
stars (Hut & Bahcall 1982). The real value for T lies somewhere 
in the middle of these extremes. 

Finally, we can calculate the distribution of relative velocities 
of single stars in the instantaneous rest frame of a binary. It is 
again a Maxwellian distribution but with an effective mass m* 
and temperature T*: 



2m 2 



T* 



m + 2m 3 

m* m* , 
— T + —T 
m 2m 



(18) 



2 1 , 
—T + -T . 

3 3 



And we have for the velocities 

,1/2 / * x 3/2 

/<>•'.> 



\kT* 



v exp - 



2kT* 



(19) 



From now on we will characterize the temperature T* by the 
thermal velocity dispersion in the relative velocities: 



=(v*n=3-. 



(20) 



Although these expressions were obtained for equal-mass bina- 
ries, the same procedure can be extended to other masses. The 
velocity dispersion is the classical parameter used to character- 
ize the relative velocities of different types of stellar systems. 
Typical values for s are extremely dependent on the environ- 
ment, and may be as low as 0.3 km/s in open clusters, or up to 
two orders of magnitude higher in the solar neighborhood (see 
Binney & Tremaine 2008 for a detailed analysis). 

From equations (26), (fT5l l, (fT9l and ( l20l i. we can calculate 
the probability function P(r m ,„, v*) dr m i„ dv* for a given value of 
s. Results are shown in Figure |6]for three magnitudes of the ve- 
locity dispersion, and adopting values for dr mm and dv* from a 
100 x 100 grid of initial conditions in the plane. 

As expected, higher velocity dispersions increase the prob- 
ability of encounters at higher values of v*. For 5=1 km/s, 
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Fig. 5. Distribution of the final inclinations of the planetary orbit 
with respect to the orbit of the binary for the full 3D simulations. 
Black histograms show all type II (left) and type III (right) out- 
comes. Those leading to planetary eccentricities below 0.1 are 
shown in red, while the accretional friendly final systems are de- 
picted in blue and green. 



only encounters with vo ~ v esc are relevant and most of the solu- 
tions discussed in the lower plot of Figure [3] will have negligible 
weight in the overall outcome probability. On the other hand, 
for .s = 10 km/s the spread in v* is much larger and even high 
incoming speeds are statistically pertinent. 

We can now incorporate these distribution functions into the 
results shown in Table [2] to give more realistic probabilities of 
positive outcomes. Before discussing these results, we must also 
incorporate a distribution function for the impactor mass. For 
this we have adopted the so-called universal initial mass function 
proposed by Kroupa (2001) which, in the mass interval between 
~ O.1M and ~ 1M , can be expressed by a multiple-part power- 
law £(m) oc nr a> with 



ai = 1.3 if 0.08 < m/M Q < 0.5 
a 2 = 2.3 if O.5O<m/M 



(21) 



and where £(m) dm is the number of single stars expected in the 
mass interval [m,m+dm]. This distribution function is shallower 
than the classical Salpeter (1955) prescription for small stellar 
masses but, at least for our application, both give similar results. 

Taking into consideration all distribution functions, we can 
now estimate the final probability for favorable scattering out- 
comes, integrated over the impactor mass mc and initial values 
for r m i„ and v*. Results are shown in Table [3] (which adds all 
positive type I + type II + type III outcomes), where we have 
considered three different possible values for the velocity dis- 
persion s of background stars. We have also performed separate 
calculations assuming that the maximum impact speeds for plan- 
etesimals leading to accretional collisions is either 10 m/s or 100 
m/s. We recall that the first is a very strict limit in which all col- 




v [v_esc] 




v [v_esc] 



s = 10 km/s 




Fig. 6. Relative probability function P(r m ,„, v*) a?r ml „ dv* for stel- 
lar encounters with given values of the incoming speed v (in 
units of v c and minimum distance r„„,„ between the impactor and 
the original binary system (in units of the apocentric distance 
Qab). In each graph we adopt a different value for the velocity 
dispersion s. 

lisions should be constructive (Thebault et al. 2006), while the 
higher impact velocity allows some disruptions to take place. 

Contrary to expectations, results are not extremely sensi- 
tive to the adopted value for s. For very stringent accretional 
restrictions, we find that around 1.5% of all stellar scattering 
events lead to favorable outcomes, while this number increases 
to ~ 5.5% if we allow higher impact speeds. This seems to imply 
that given a "normal" planetary system consisting in either (i) a 
planet around a single star interacting with a binary stellar sys- 
tem, or (ii) a planet in a circumstellar orbit around a wide binary 
interacting with a single star, the final result could be a system 
similar to y-Cephei in (roughly) 1% - 5% of cases. 



6. Conclusions 

The discovery of at least four exoplanets in circumstellar orbits 
in tight binary systems has raised the question of their formation 
history. So far, all attempts to explain their origin via normal in- 
situ formation processes have been unsuccessful, and we won- 
der whether other more exotic mechanisms may have played a 
role. Following the suggestion of Marzari & Barbieri (2007a), in 
this work we have analyzed whether these systems, in particular 
y - Chepei, could have begun their lives in other dynamical con- 
figurations more suitable for planetary formation, attaining their 
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Table 3. Percentage of positive outcomes for two limiting im- 
pact velocities Av for accretion and three different velocity dis- 
persions s for background star. Values are integrated over dis- 
tribution probabilities for stellar mass, impact speeds and cross- 
sections. 



present structure later on as a consequence of a close encounter 
with a background singleton. 

Following a series of numerical simulations of time-reversed 
stellar encounters, whose results were later weighted adopt- 
ing classical expressions for probability distributions for stel- 
lar mass, relative velocities and impact cross-sections, we have 
found that between 1 and 5 percent of fly-bys between a plane- 
tary system and background stars could lead to planets in tight 
binary systems. Although this number may not seem high, it 
could in fact explain why planets in tight binaries are not more 
numerous, which would indeed be the case if in-situ planetary 
formation were possible. 

However, we must remain cautious. Although we have at- 
tempted to give some estimate about the probability outcomes, 
given a sufficiently large number of initial conditions and val- 
ues in the parameter space, it is possible to obtain almost any 
result from a scattering event. This does not mean that the event 
actually happened or that the explanation lies there. However, 
as Sir Arthur Conan Doyle said through his character Sherlock 
Holmes: "Once you eliminate the impossible, whatever remains, 
no matter how improbable, must be the truth." 
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